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Abstract. Shock waves associated with shell type supernova 
remnants are considered to be possible sites of cosmic ray ac- 
celeration. Since shocks are capable of accelerating electrons 
in addition to protons one anticipates both species to con- 
tribute to the high energy radiation expected from these ob- 
jects. Adopting a simple model for particle acceleration we 
calculate in a self-consistent manner the time-dependent syn- 
chrotron and inverse Compton radiation of high energy elec- 
trons assumed either to be accelerated directly by the shock 
wave or to be injected at high energies as secondaries from the 
hadronic collisions of relativistic protons with the circumstellar 
material. We deduce that for standard supernova parameters 
the TeV flux produced from neutral pion decay is about the 
same order as the flux expected from directly accelerated elec- 
trons. 

Key words: acceleration of particles - Radiation mechanisms: 
Compton and inverse Compton - Radiation mechanisms: cy- 
clotron and synchrotron - shock waves - cosmic rays - super- 
novae: SN 1006 



1. Introduction 

Shell-type supernova remnants (SNRs) are long thought to be 
the sources of the nuclear component of cosmic rays up t o en- 
ergies close to the k nee ( see, for example, Axford [1981] and 
Lagage & Cesarsky [1982]). SNRs appear to be one of the few 
galactic objects which release enough energy to satisfy the ob- 
served flux of cosmic rays. Also recent developments in the 
theory of diffusive acceleration in shock waves provide the re- 
quired theoretical b ackgro und for such a picture (see the review 
by Jones & Ellison Il99l|). 



In addition to protons, electrons are also expected to be ac- 
celerated directly in these shock waves. However, as protons are 
expected to carry practically all the energy available for accel- 
eration (usually taken to be 10% of the total supernova energy) 
the energy content in accelerated electrons is not known. Nev- 
ertheless, if one assumes that the measured ratio of cosmic ray 
electron energy density to that of protons is the same as the 
one which is produced by the SNRs then we should expect that 
about 1% of the SNR energy available for acceleration goes to 
electrons. 



Despite the theoretical arguments, there is not as yet any 
firm observational evidence to support the above picture of 
cosmic ray origin. It was suggested that a direct proof for 
the theory of proton accelera tion w ould be the detec tion of 
SNRs at TeV energies (Dorfi ]l99l|] and Drury et al. pjgl- 



henceforth DAV94). Relativistic protons undergo hadronic col- 
lisions with the ambient matter producing TeV 7— rays from 
neutral pion decay. This flux should be, in principle, detectable 
with present-day Cerenkov detectors. 

On the other hand, the re is s trong observati onal e vidence 
that SN 1006 (Becker et al. |l98f|, Koyama et al. |l995|]) shows 
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a non-thermal X-ray spectrum and the possibility that the 
accelerated (primary) electrons can radiate synchrotr on X- rays 
has already been discu ssed by Reynolds & Chevalier ( 1981 ) and 
Ammosov et al. ( 1994 ) -however see Hamilton et al. ( |1986| ) for 
an alternative explanation. It has been suggested that if the X- 
rays observed from SN 1006 prove to be non-thermal in origin, 
this can only mean that electrons are accelerated to energies 
of 100 TeV or more, providing thus for the first time evidence 
for very high energy particle acceleration in SNRs. 

The above statements, i.e. that TeV radiation is proof of 
proton acceleration while the non-thermal X-rays are an indi- 
cation for electron acceleration, ignore the fact that each of the 
accelerated species can give independently a contribution to 
the X-ray and 7— ray flux. Thus directly accelerated electrons 
will not only produce synchrotron X-rays but also 7— rays 
from inverse Compton scattering on any ambient photon field 
present. Similarly, accelerated protons will not only produce 
neutral pions but also charged pions which, upon decay, will 
create secondary electrons. These electrons will radiate just as 
any directly accelerated electron contributing to the total SNR 
X-ray and 7— ray flux. The question which arises therefore is 
the following: If SNRs turn out to be sources of TeV radiation 
will the radiation come mostly from protons (ir° decay) or from 
electrons (inverse Compton scattering)? And if there is, after 
all, a non-thermal component in the X-ray regime will it be 
mostly due to directly accelerated electrons or due to secondary 
electrons produced in hadronic collisions? The above questions 
are obviously of some importance to acceleration theories. 

The aim of the present Letter is to try to address the ques- 
tions posed above. In doing so we calculate, within the stan- 
dard SNR particle acceleration framework, the time-dependent 
X-ray and 7— ray fluxes produced by each of the accelerated 
species. In §2 we describe the model we base our calculations 
on, in §3 we give results and we conclude in §4. 
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Fig. 1. Maximum energy of electrons as a function of the magnetic 
field at shock B. £ a cc is the maximum energy the particles can attain 
when losses are neglected - see text. Ei oss gives the critical energy 
at which the acceleration rate can be balanced by the energy loss 
rate. Thus the maximum possible energy is given by the solid line. 
Ei and Eio are the energy electrons should have in order to radiate 
synchrotron X— rays of energies 1 and 10 keV respectively 

2. The Model 

We base our SNR model on a general picture including diffu- 
sive a cceleration (for detail s of the underlying physics see Dorfi 
[1993 1 and Kirk et al. [1994]) treating both electons and protons 
as test particles. We assume that the SNR shock front is moving 
through the interstellar medium with a velocity ~ const, in 
the free expansion phase and u s h oc t~ 3 ^ 5 in the Sedov phase. 
The Sedov phase starts approximately at a time tsd when the 
shock has reached a radius such that the swept up mass equals 
the ejected mass, thus tsd — 2W£l^ > sl pi 1 '' 3 u~ h '{ t 3 yr. Here 
£sn,51 is the total energy released in the supernova explosion 
in units of 10 51 erg, pi is the external matter density in units 
of 1 H-atom/cm 3 and « s h,4 is the shock velocity in units of 10 4 
km/sec. 

We take a simple Bohm type of diffusion upstream and 
downstream of the shock and we assume that protons as well 
as electrons are accelerated there. We determine the max- 
imum energy to which the particles have been accelerated 
after some time t by simply integrating over the accelera- 
tion rate (see, for example, Webb et al. [ 1984 J) which gives 



£ acc ~ 3A(r c )B_5U sh 4 t yr TeV for t < t Sd and _E acc ~ 
630yl(r c )B_5U s 1 / 3 4 £: s 1 ^ 3 51 p 1 " 1/3 TeV ~ const for t > t sd . Here 
A(i c ) = (r c — l)/r c (r c + 1) with r c being the compression ratio 
of the shock, while _B_s is the magnetic field strength given in 
units of 10 -5 G. Especially for the case of the accelerated elec- 
trons the maximum energy the particles obtain can be limited 
by losses. We estimate therefore the critical energy Ei OBB (t) 
at which the acceleration rate can be balanced by the loss 
rate (synchrotron and inverse Compton combined) and set 
Ee,nmx(t) = min[i5 acc (t) , E loss (i)] and this is shown in Fig. 1 
as a function of the magnetic field B. 

2.1. Emission due to primary accelerated electrons 

We calculate the radiated spectra from the relativistic elec- 
trons once these have escaped downstream of the shock. We 



treat this problem in a self-consistent manner by solving first 
a time dependent kinetic equation for the relativistic electron 
distribution function n e (E e ,t) 



dn c (E c ,t) 
dt 



Q e pm (E e ,t)+C e (E e ,t). 



(1) 



Here €. a (E e , t) denotes the electron losses which are taken to be 
due to synchrotron radiation and inverse Compton scattering, 
as well as losses due to adiabatic expansi on u ( for the relevant 
expressions used see Mastichiadis & Kirk |199q ] and references 
therein). The source of high energy electrons is given by the 
rate Q prm of primary electrons that escape downstream. For 
them we assume a power law spectrum in energy and write 
Q| rm (-E e ,t) = Q efi (t)E~ s with E c < E e , max (t) and s = (r c + 
2)/(r c — 1). This approximation is known as the "Onion-shell- 
model" (Bog dan fc Volk 1983). For determining Q e .o(t) we 
follow Drury and assume that a part of the total energy 

flux through the shock goes into accelerating electrons. Thus 
we write 



— R(t) 3 I dE c E c Q c , pTln (E c ,t) 



^[R(t)Y- P [u sh (t)Y 



(2) 



where p is the ambient matter density and £ is a proportional- 
ity constant to be determined. For steep electron spectra (i.e. 
s > 2) as the ones we will examine here this energy flux is oc t 2 
for t < tsd and oc t^ 1 for t > t sd . From this form of source 
function it becomes evident that the bulk of energy output 
occurs at times t ~ fed- 

The solution of the electron kinetic equation (Eqn. [lj) gives 
the distribution function of primary electrons downstream 
n e (E e ,t) which is an implicit function of the quantity £ (Eqn. 
|2|). We normalize n e (E e ,t), and thus determine £, to the frac- 
tion of the total supernova energy that goes to accelerated 
electrons ?7 e i through the equation 



yiSsNR / dE e E e ne{E e ,tsisiR) = £ c i = ??ei£sN. 



(3) 



Here £snr is the time at which the remnant ends its life, 
7?snr = R(tsNR), while £ e \ is the energy content in elec- 
trons produced in the lifetime of the supernova. Since most 
of the power in particles is put into the SNR at times tsNR, 
at late epochs Eqn. (^) changes very slowly with time and 
therefore it is insensitive to the exac t valu e of £snr. How- 
ever for concreteness we follow Dorfi ( 1993| ) and use tsNR 



1.3 lO^^gJ/^pJ' 1 ' 5 '' 33 y ears - Once n e (E e ,t) is determined we 
can use it to calculate both the synchrotron and inverse Comp- 
ton emissivities by folding it with the corresponding emissivi- 
ties. 

2.2. Emission due to accelerated protons 

In the case of proton acceleration we solve a kinetic equation 
similar to Eqn. (|l]) but taking into account only adiabatic 
losses. The source of high energy protons is provided by the 
proton escape downstream. In close analogy with the electron 
case we introduce a proton efficiency r) pr which we relate to the 



In dense media electron bremsstrahlung will become impor- 
tant, however we restrict our analysis here in shocks propagat- 
ing in the ISM where the densities involved allow us to neglect 
the above process. 
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Fig. 2. Synchrotron spectra produced due to primary accelerated 
electrons at t= 1,000 years after the explosion as a function of the 
magnetic field, (a) is for 3 /tG, (b) is for 10 /tG, (c) is for 30 fiG and 
(d) is for 100 /iG. The other parameters are as given in the text. 



total energy content in protons £ pl produced in the lifetime of 
the supernova by an equation similar to Eqn. The solu- 
tion of the kinetic equation for protons gives their distribution 
function n p (E p ,t) which we use subsequently to calculate the 
produced spectra of the 7— rays Q wl and secondary electrons 
Ql ec (E c ,t) (produced through neutral and charged pion decay 
respecti vely - for details on the methods used see Mastichiadis 
fc Kir k [1995]). Q,r 7 is the 7-ray emissivity considered by Dorfi 
( tl99l| ) and DAV94. 

For the case of secondary electrons we have to solve again 
the electron kinetic equation ([j]) using Qs C c(E c , t) instead of 
Q p[m {E e ,t) as the source of energetic electrons to obtain the 
distribution function of the secondary electrons n e]aec (E e ,t). 
This can be used subsequently to calculate both the syn- 
chrotron and inverse Compton emissivities in a manner similar 
to the one described in Section 2.1. Note that once n p (E p ,t) 
has been specified, the quantities Q„- 7 , Qg CC and n GiSec are com- 
pletely determined. 
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3. Results 

We solve numerically the above described three kinetic equa- 
tions (two for electrons - primary and secondary, and one 
for protons). We assume that the shock compression ratio is 
r c = 3.73 which implies that the accelerated particle index is 
s — 2.1. This value was chosen because it is in agreement with 
estimates for the possible co smic ray index at source - see 
DAV94 and Berezhko et al. (1994). Furthermore we assume 



that r] pr — .1 and r) e \ = .001 while Esn = 10 51 erg, u a h = 10 4 
km/sec and p=l H-atom/cm 3 , while we leave the magnetic 
field strength as a free parameter. Note that our results scale 
linearly with r] pl and r] e i. Also the proton related fluxes (i.e. 
Q„- 7 and Qs CC ) scale linearly with the ambient density p. 

Fig. 2 shows the synchrotron X-ray spectra of primary elec- 
trons at t=l,000 years after the explosion as a function of the 
magnetic field which we assume to be uniform in the down- 
stream region. One notices that both the flux and the high 
energy cutoff increase with increasing magnetic field. The flux 
increases because of the dependence of synchrotron emissivity 
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Fig. 3. A SNR synchrotron X— ray lightcurve in the 1-10 keV energy 
band. The full and dashed lines are due to radiation of primary elec- 
trons in a magnetic field of B=100 /iG and 10 /tG respectively, while 
the dotted line is due to secondary electrons and B=100 /tG. B SNR 
7— ray lightcurve for energies > 1 TeV. The full line is due to 7T° de- 
cay, while the dashed and dotted lines are due to inverse Compton 
scattering of primary electrons on the microwave background. The 
magnetic field was taken to be 3 and 100 /iG respectively. 



on the magnetic field. The high energy cutoff increases because 
one can roughly write e ma x oc E^ m!ix B oc B 3 for those magnetic 
fields where acceleration is not limited as yet by losses (see Fig. 
1). From Fig. 2 it is apparent that one needs to have magnetic 
fields larger than 3 /iG for 10 keV X-rays to be produced (see 
also Fig. 1). There are some interesting consequences to the 
spectra for relatively high magnetic fields (B> 30 /iG). Thus 
the high energy cutoff stops increasing as now -E e ,max is deter- 
mined by losses. Also the photon spectrum shows a break; this 
reflects a similar break in the electron distribution function 
which is caused by electron synchrotron cooling. This break is 
absent for low magnetic fields as for these fields the synchrotron 
timescale is long compared to the age of the SNR. 

Fig. 3 shows the expected X-ray and 7-ray (> 1 TeV) 
lightcurves from a SNR characterised by the parameters given 
above. X-rays are produced by the synchrotron radiation of 
electrons while 7-rays are produced from n° decay and from 
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electron inverse Compton scattering on the cosmic microwave 
background. For the assumed values of n c \ and n pl we find that 
primary electrons produce higher fluxes than secondaries both 
in the X- and in the 7— ray regime. For high magnetic fields ( > 
50 fiG) the synchrotron cooling timescale is short compared to 
the age of the SNR; this means that electrons radiate efficiently 
and the synchrotron lightcurve follows essentially the source 
function of the electrons (solid line curve in Fig. 3A). When 
the magnetic field is lower the above argument is no longer 
valid as the electrons never have the time to cool efficiently 
(dashed line curve in Fig. 3A). 

In the 7-ray regime (Fig. 3B) we find that the inverse 
Compton contribution from the primary electrons (dashed line) 
is of the same order as that of the it 7-rays (solid line) pro- 
vided that the magnetic field is not much larger than 10 /iG. (It 
is interesting to note that our tv° 7-ray lightcurve is in excellent 
agreement with the result of DAV94 despite the different ap- 
proaches of the two papers.) Higher magnetic fields suppress 
the inverse Compton flux as they put most of the radiated 
power in synchrotron. 

However we find that TeV observations could tell whether 
the observed spectrum is due to protons or electrons. In the 
case of protons one expects a spectrum with a spectral index 
of s simply because the produced pions reflect the spectrum 
of the relativistic protons. In the case of electrons, however, 
the spectrum should be Hatter having an index a — (s + l)/2 
-assuming that we are away from the cutoff. This happens 
because the electron distribution function should not be steep- 
ened by losses and therefore will have a slope of index s which, 
in turn, gives the aforementioned index in radiation. Thus tak- 
ing the usual values for s between 2.1 and 2.3, one finds a to 
vary between 1.55 and 1.65. 

Finally, it is interesting to note that, according to the 
present model, the >1 TeV flux from SN 1006 due to primary 
electrons will be ~ 2 10 -12 ph/cm 2 /sec which makes this ob- 
ject detectable with present day Cerenkov detectors. At the 
same time the low number d ensity around SN 1006 (.05 H- 
atoms/cm 3 - Hamilton et al. [1986]) makes the corresponding 
flux due to 7T° 7— rays about an order of magnitude less. 



4. Conclusions 

Adopting a simple model for diffusive acceleration at shock 
fronts we have calculated the non-thermal luminosity of SNRs 
by treating the radiation problem in a self-consistent manner; 
thus we first solved the time dependent kinetic equations for 
electrons and protons and consequently we used the obtained 
distribution functions to calculate the radiated spectra. We 
have shown that even for ambient matter densities as high as 
1 H-atom cm -3 , the X-ray luminosity produced by the syn- 
chrotron radiation of secondary electrons is many orders of 
magnitude below the radiation of the directly accelerated (pri- 
mary) electrons provided that the latter are given about .001 
of the total SN energy budget. Thus a possible confirmation of 
the non-thermal char acter of the X-ray flux detected from SN 
1006 (Koyama et al. [ 1995 1 ) will indeed indicate that the radi- 
ating electrons are directly accelerated by the SNR blast wave 
and are not injected as secondaries from hadronic collisions. 

On the other hand we find that the directly accelerated 
electrons can produce a TeV flux by inverse Compton scatter- 
ing on the microwave background which is about of the same 
order as the flux produced from the n° decay (DAV94) pro- 



vided that the magnetic field is not much higher than 10 ^tG. 
Therefore a possible detection of a SNR at TeV energies will 
not neccessarily mean that protons are accelerated there. To 
deduce that one needs detailed spectral information as elec- 
trons are expected to give flatter power laws. 

Note also that as we have assumed that the high energy 
electrons scatter only on the microwave background, all the in- 
verse Compton fluxes presented here are only lower limits since 
the inclusion of any other photon field ( such a s the galactic in- 
frared background-see Cox & Mezger [ 1989 1 ) would enhance 
the inverse Compton emissivity. Interestingly enough, TeV ob- 
servations of SNRs have started co nstra ining the theories of 
proton acceleration (Lessard et al. [1995]), so this constraint, 
when applied to electron radiation alone, might be used to ob- 
tain a limit on the magnetic field (De Jager et al. [ 1995 ]). 

Finally, one should add that our results do not simply scale 
with the ambient matter density when this becomes large. In 
such a case electron bremsstrahlung becomes important not 
only as a radiation mechanism but as an energy loss process as 
well. This effectively is an additional term in Eqn. (Q) which 
will affect the electron distribution function and thus the radi- 
ated spectrum. 
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